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We show that a two-dimensional spin-orbit-coupled system in the presence of a charge/spin-density wave with 
a wave- vector perpendicular to an applied electric field supports bulk manifestations of the direct/inverse spin- 
Hall effect. We develop a theory of this phenomenon in the framework of the spin diffusion equation formalism 
and show that, due to the inhomogeneity created by a spin-grating, an anomalous bulk charge-density wave is 
induced away from sample boundaries. The optimal conditions for the observation of the effect are determined. 
The main experimental manifestation of the bulk spin-Hall effect, the induced charge/spin-density-wave, is 
characterized by a ;r/2-phase shift relative to the initial non-homogeneous spin/charge-polarization profile and 
has a non-monotonic time-varying amplitude. 



Electron spin transport in semiconductors has recently 
evolved into a subject of intense research as key element of the 
rapidly developing field of spintronics (Tj2]. One of the main 
challenges is to generate spin polarization and to transport 
spin in nonmagnetic materials using electric fields, by taking 
advantage of the coupling between spin and orbital degrees of 
freedom |3 , 4|. Of particular importance in this respect are a 
family of anomalous transverse transport phenomena, such as 
the spin Hall effect I? '51, which has received recently tremen- 
dous attention. The experimental manifestation of the effect 
is the appearance of spin accumulation near the edges of the 
sample|l2l|8]|9l, if an electric current is driven through a sys- 
tem with either intiinsic or extrinsic (i.e., impurity-induced) 
spin-orbit coupling flOl . Hence, the canonical spin-Hall ef- 
fect is, at least from an experimental perspective, an edge phe- 
nomenon whose magnitude depends strongly on the specific 
boundary conditions lflTlfT2l . The role of the edge is to create 
a strong inhomogeneity where the experimentally observable 
spin density can accumulate. However, it is possible to cre- 
ate large length scale inhomogeneities in a controlled way, for 
example by generating a modulated charge or spin distribu- 
tion or a spatially varying spin-orbit coupling. By analogy 
with the canonical edge spin-Hall effect, the externally gen- 
erated charge/spin-densities would effectively create multiple 
boundaries in the bulk, thus opening the possibility of realiz- 
ing bulk manifestations of the direct/inverse spin-Hall effect 
and allowing for a controlled study of this phenomenon. 

An effective way of producing bulk manifestations of the 
spin-Hall effect, suggested by the recent work of Koralek 
et al. \\13l . is to use the transient spin grating (TSG) tech- 
nique lfT4l [TSl [T6l to generate and monitor time dependent 
spin and charge profiles. Within the TSG method, a sinusoidal 
spin-polarization wave is generated by two interfering non- 
collinear laser beams with orthogonal linear polarization. This 
induces a modulation in the index of refraction, which can be 
measured at subsequent times by the diffraction of a probe 
pulse. In the presence of an external electric field oriented 
perpendicular to the spin-polarization wave-vector, a charge 
density with the same wave-vector is expected to develop (see 
Fig. [TJ. Alternatively, in a spin-orbit interacting system a 
charge density wave is expected to induce spin modulation 
under the action of an external electric field. 



In this Letter, we develop a theory of the bulk spin-Hall ef- 
fect in the diffusion Umit, in the presence of Rashba H and 
(linear and cubic) Dresselhaus f3l spin-orbit interactions. We 
focus on the time evolution of a charge density profile in- 
duced by an optically generated spin-polarization wave and 
its dependence on the spin-orbit couplings and on the spin- 
grating wave-vector. In particular, we determine the optimal 
parameters for observing the spin-Hall effect with spin grat- 
ings. These optimal parameters result from a balance between 
two competing requirements: 1) To create slowly decaying 
spin-polarization waves, and 2) To have a strong spin-charge 
coupling. The first requirement is related to the more general 
challenge in the field of spintronics of identifying mechanisms 
allowing for long spin relaxation times. In the presence of dis- 
order, spin-orbit interactions lead to spin relaxation through 
the Dyakonov-Perel mechanism [17 ,1 . Recently, it was shown 
that an enhanced spin life time can be realized by tuning the 
spin-orbit coupling so that the Rashba and the linear Dressel- 
haus couplings become equal |fT3l. In this regime, SU(2) spin 
symmetry is restored, allowing for a long lifetime helical spin 
density mode termed the "persistent spin helix"| 18|, provided 
that the cubic Dresselhaus contribution be minimized |19|. 
However, in the persistent spin helix regime the coupling be- 
tween the spin and the charge channels vanishes and the spin 
Hall effect cannot be observed. Hence, the second require- 
ment: the existence of a strong spin-charge coupling. 

We consider a two-dimensional electron gas in a III-V type 
semiconductor quantum well grown along the [001] axis (set 
as the z axis). In the presence of Rashba 111, as well as linear 
and cubic Dresselhaus (3] spin-orbit interactions, the Hamil- 
tonian describing the conduction band electrons is 

>H= f-+h(p).o- (1) 

2m 

where m is the effective mass, & = (o"a, o"i/) are Pauli ma- 
trices and h(p) - {h^,hy) is the momentum-dependent ef- 
fective "magnetic" field describing the spin-orbit interaction. 
The Rashba contribution, h*(p) = avf{-py, p^), with Uf the 
Fermi velocity, arises from the inversion asymmetry of the 
quantum well confining potential, with the coupling constant 
a measuring the strength of the Rashba spin-orbit coupling 
relative to the Fermi energy. In addition, the lack of inversion 
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FIG. 1: (Color online) Charge-density profile induced by the relax- 
ation of a spin-density wave in the presence of a uniform electric 
field. The initial spin density corresponds to a sinusoidal wave with 
wave- vector q of the out-of-plane 5; component, as symbolized by 
the blue (spin up) and red (spin-down) arrows. The in-plane electric 
field is oriented perpendicular to q. Notice the f shift of the induced 
charge-density profile relative to the spin-density wave. 

symmetry of the semiconductor crystal gives rise to the Dres- 
selhaus spin-orbit interaction, h^'(p) = PyVpiPx^-Py) ™d 
h^'(p) = -4(]2^-(Pxpl,-PijP])y where pf is the Fermi mo- 
mentum and fix and fij, are dimensionless coupling constants 
for the linear and cubic Dresselhaus terms, respectively. In the 
presence of disorder, the coupled spin and charge dynamics 
can be described by a generalized diffusion equation, which in 
the absence of an external electric field has the form liT9ll20ll 

{d, - 2>V^)pi = (r'^' - P'^*Vt + ■ V)pj, (2) 

where po is the charge density and p 1,2.3 = Px.y.z ^re spin 
densities. The parameters F'-i describe the Dyakanov-Perel 
spin relaxation IfTTI . ^ = is the diffusion constant, 

with T the mean scattering time, P''* = -p>'* characterize 
the precession of the spin polarization and and C'-' describe 
the coupling between the spin and charge degrees of free- 
dom. In momentum space, the diffusion equation becomes 
p; = 0, where n,y have coefficients given by 
pj, P'}^ and C'^ US). The formal solution of the diffusion 
equation is p,(r, f) — J dr'Dij(r,r',t)pj(r',0), where p,(r, 0) 
is the initial spin-charge distribution and D = [1 - O] ' is the 
Green's function of the diffusion equation, or the diffuson. 

The generalization of the spin-charge diffusion formalism 
developed in Ref. 19 for the case of a uniform electric field 
amounts to the formal substitution |21 1 

V ^ V + (3) 

where £ is a uniform electric field and yu is the mobility of the 
two-dimensional electron gas. Note that, neglecting the spin- 
charge coupling, this substitution generates the standard drift- 
diffusion equation for the charge channel, while the descrip- 
tion of the spin sector is in agreement with a semi-classical 
kinetic theory of electron spin transport derived using the 
Keldysh Green's function formalism [21] . The substitution 



^ is valid as long as non-linear contributions of order O(E^) 
are small and assuming that the effects of the electron-electron 
Coulomb interaction can be neglected. Without loss of gen- 
erality, we focus on the geometry corresponding to Fig. [T| 
and consider a system with an initial out of plane spin-density 
wave, pj(r, 0) - n()Cos(qr+), oriented along the [110] direction 
(e+) and a weak constant electric field, E = EqC-, oriented 
along [110] (e_). In momentum space, the substitution ^ 
becomes k — > k - i/uE/^^ and the inverse of the diffuson is 
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where E - iiELJlS! is a dimensionless measure of the elec- 
tric field strength and s - -iio{q) + q^ + \. All lengths are 
measured in units of spin relaxation length, L, - l/2pfA and 
times in units of spin relaxation time, t, - Ir/g'^A, where 

A - (a^ + (J3i - Pi)^ + pfj ^nd g - IvfPft is a dimension- 
less conductance. The spin-spin coupling parameters [19], 
T± = Tr ± 7d^ with y« = 2a/A and yo - 2(Pi -y63)/A, are in- 
dependent of the overall strength of the spin-orbit interaction 

r = (^a^ +y6i + ySj) ^ and lie within a disc of radius 2. The 
spin-charge coupling parameters 1 19J, A± = Ai +A2, with Ai = 
[(3;63-y8i)(ff'-y6?+j6^)-j6iy8^]/Aand^2 = a(a^-/3l+6/3l)/A, 
are quadratic in the spin-orbit interaction strength. 

The induced charge density po(r, t) is determined by the 
matrix element D03 = Iiio^'(<J)*^"^''^'''' °f *s diffuson. Here 
icoiiq) are the relaxation modes obtained from the equation 
det |i - fl(w, q)j - and A/(q) are momentum-dependent 
amplitudes. To order 0{E^) the relaxation times are indepen- 
dent of the electric field, while the amplitudes have a linear de- 
pendence, A/(q) - i(fj.E/vf)qAi(q), with Ai{q) being an even 
function of momentum and q = q-e+. Note that, if one initially 
generates a charge density profile, the external electric field 
induces a spin wave with a spatial and time dependence de- 
termined by Z)3() - -Z)o3. Hence the present analysis applies 
to both the direct and the inverse spin Hall effect. Explicitly, 
an initial out-of-plane spin density wave pj(r, 0) = nocos(^r+) 
induces a time-dependent charge density wave 

p(r, f) = no sin(^r+)^ V qAiiq)e-''^'^'^")' (5) 

Note that the induced charge density wave (CDW) is phase 
shifted by :7r/2 relative to the initial spin-density wave (see 
Fig. [T]i and has a time dependent amplitude noifiE /uf)A{t) 
where A(r) = qAi{q)e''"'^i'^' . The general behavior of 
the induced CDW amplitude A(f) is shown in Fig. |2] At f = 
the amplitude of the CDW vanishes, as the system is initially 
uniform, while at long times A{t) decays exponentially with 
a characteristic lifetime \ liiuiiq)) given by the lowest fre- 
quency relaxation mode. At intermediate times of order the 
CDW amplitude has one maximum and/or one minimum. The 
largest absolute value defines the peak amplitude A„„„. The 



FIG. 2: (Color online) Time dependence of the induced charge- 
density wave amplitude A(t) for various values of the dimensionless 
spin orbit coupling parameters {yR,yD) and for an overall spin-orbit 
coupling strength T = .001. The wave-vector q || e+ has a fixed 
value, q = 0.6/L,. The amplitude of the induced spin-density wave 
varies non-monotonically and is characterized by a peak value A„„,. 
and an exponential decay at large times. 

strength of the Rashba and Dresselhaus spin-orbit interaction 
in GaAs quantum wells can be adjusted by varying the doping 
asymmetry or the width of the quantum wells. Values in the 
range of Of = 0.5 x 10"^ H- 1.5x10"^ and/?; = lxlO"^H-3xlO"^ 
withySs = 0.3x10-^ till can be experimentally achieved, thus 
most of region in the vicinity of the boundary of the radius 2 
disc in the (y^, jo) parameter space can be probed. Scaling a, 
/3i and ^63 equally will not change the spin-spin couplings ys 
or yj), but it will change the spin-charge couplings A+ and /l_ 
which are quadratic in the overall spin-orbit coupling strength 
r. The amplitudes Ai(q) depend linearly on A+ and /i_ with 
higher order corrections of order A^.. Thus for experimentally 
realizable two-dimensional spin-orbit interacting electron sys- 
tems characterized by F <K 1, the higher order corrections 
due to the spin-charge couplings are negligible and the ampli- 
tude A(f) is approximately linear in the spin-charge couplings. 
Since the factor of q in A(f) gives a contribution of l/F, as 
the wave-vector is measured in units of 1/Lj, we conclude 
that the amplitude A(f) of the induced charge-density wave 
depends linearly on the overall spin-orbit interaction strength 
r. This proportionality relation holds as long as we express 
the wave-vector in units of I/L.5. Furthermore, we find that 
the induced CDW amplitude A(f) is independent of the dimen- 
sionless conductance g, provided time is measured in units of 
l/r,. Consequently, the bulk manifestation of the spin-Hall 
effect proposed here, can be enhanced by reducing the carrier 
density of system, which will increase the ratio between the 
strength of the spin-orbit interaction and the Fermi energy. 

Next, we study the dependence of the induced charge den- 
sity wave amplitude on the ratio between various components 
of the spin-orbit interaction for a fixed value of the overall 
spin-orbit coupling strength F. Fig. |3] shows the maximum 
amplitude of the charge-density wave, A,„„^, for the experi- 
mentally relevant spin-orbit coupling strength F - 0.001 and 
wave-vector q = 0.7/Lj. The peak amplitude vanishes for 



FIG. 3: (Color online) Dependence of the absolute value of the 
peak amplitude A„,„j- on the spin-orbit parameters (70, y^) for a fixed 
value of the overall spin-orbit interaction strength, F = .001, and 
q = 0.7/L,. A,„aj. vanishes for pure Dresselhaus spin-orbit cou- 
pling, 7r = (the segment between the horizontal pair of white 
dots), pure Rashba coupling, (y^ = 0,7^. = ±2) (vertical pair of 
white dots), and at the symmetry points (70 = ± V2, = ± y/l) 
(green dots). The maximum of the peak amplitude corresponds to 
(70, 7;;) = (-1.08,-1.25) (inside the lower left quarter of the param- 
eter space, A,„a j = 7.8 x 10"*), while three other local maxima are lo- 
cated at (y^, 7^) = (-1, 1.24) (upper left quarter, A„,„v = -5.5x10"''), 
(0.80, 1.06) (upper right, A„„ = 3.2x 10"''), and (0.76, -0.98) (lower 
right, A,„„j = -2.8 x 10"''). All these maxima involve large relative 
contributions of the cubic Dresselhaus coupling, ySj /F = 0.5 -r 0.68. 

pure Dresselhaus spin-orbit coupling, y^ = 0, pure Rashba 
coupling, (yo - 0, y^ = ±2), and at the symmetry points 
(yo - ± V2, yn - + V2) which support the persistent spin he- 
lix modes (see Fig. [3]). This is consistent with previous results 
showing that, at least in uniform and stationary conditions, the 
spin Hall conductivity in systems with pure Rashba or pure 
linear Dresselhaus spin-orbit interaction vanishes |22|. Our 
analysis reveals the absence of any manifestation of the spin- 
Hall effect for these types of spin-orbit interactions in non- 
uniform systems and under time-dependent conditions. The 
absolute maximum of the peak amplitude, A,„a_^ = 7.8 x 10"'*, 
is realized for {yo^yR) - (-1.08,-1.25). The corresponding 
original spin-orbit couplings are (a,/3i,/33) - (-7.4, 0.3, 6.7)x 
10"'*. Several other local maxima (minima) can be identified 
throughout the parameter space (see Fig. [3]). To enhance the 
peak amplitude of the induced charge profile, one has to con- 
sider systems with strong cubic Dresselhaus spin-orbit cou- 
pling. This condition is opposite to that required for the real- 
ization of the persistent spin helix mode lfT3l 11911 . Note that 
the diagram in Fig. [3] has no particular symmetry, as a re- 
sult of the nontiivial dependence of the spin-charge coupling 
parameters A±- on the spin-orbit couplings. 

We consider now the case of a fixed cubic Dresselhaus cou- 
pling in the range /Sj, =2x1 0"** H-4 x 1 0"'*, which is experimen- 
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FIG. 4: (Color online) Dependence of the peak amplitude on the 
linear Dresselhaus spin-orbit coupling for various values of the cubic 
Dresselhaus and Rashba couplings. The arrows mark the values of 
/3i where the Amaj changes from an absolute minimum of A(t) to an 
absolute maximum (see also Fig. [2](. 

tally relevant for GaAs quantum wells. The dependence of the 
peak amplitude on the tunable parameters a and /5i is shown 
in Fig. |4] We stress that both the absolute value and the sign 
of the spin-orbit coupling constants are important in determin- 
ing the strength of the spin-Hall effect. Finally, we note that 
the peak amplitude also depends on the wave-vector q. A,„ax 
vanishes in the limits ^ — » and q ^ oo and is maximized in 
the range 0.5 < qLs < 0.7. Increasing the spin-orbit interac- 
tion strength enhances the bulk spin-Hall effect, provided it is 
observed at larger wave-vector values. 

For completeness we note that, if the initial spin-density 
waves have an arbitrary orientation of the q-vector, a charge- 
density wave is induced even in the absence of an external 
electric field. However, this wave is in-phase with the initial 
spin wave. Adding an external electric field perpendicular to 
the wave-vector induces an additional charge density compo- 
nent characterized by a 7r/2 phase shift, as described above, 
and causes the spin and charge profiles to drift along a direc- 
tion parallel to the q-vector, i.e., perpendicular to the electric 
field. The induced charge-density wave has the form 



3 

p(r,f) = noJ]e-"''^'i^'l^a,(q)cos^q-r + {qxE)^Cl,{q)t 



1=0 



+ (q X e) A/(q) sin q ■ r -h (q x £)_ Q/(q)f 



, (6) 



where fl/(q) are the amplitudes of the in-phase charge compo- 
nent and w/(q) are the corresponding frequencies. The elec- 
tric field induces out-of -phase waves with amplitudes A;(q) 
and generates oscillatory components of the relaxation modes 
proportional to Q/(q). 

In summary, we show that a non-homogeneous spin-orbit 
interacting system supports bulk manifestations of the spin- 
Hall effect. We extend the spin-charge diffusion equations to 
the case of a constant electric field and use this tool to char- 
acterize the charge density wave induced by an initial spin 



density wave that relaxes in the presence of an external elec- 
tric field perpendicular to the spin-polarization wave-vector 
We find that the induced charge profile is characterized by the 
same wave-vector as the spin-density wave but has a phase 
shift of +7t/2. The amplitude of the induced charge-density 
wave varies non-monotonically in time and is characterized by 
a peak value and an exponential decay at large times. We show 
how to maximize the effect by tuning the relative strengths of 
the spin-orbit interactions. Finally, we mention that similar 
non-homogeneous perturbations may lead to bulk manifesta- 
tions of the topological quantum spin-Hall effect Ii23t i24J in 
spin-orbit interacting insulators |25|. 
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